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THE  LINEARIZED  CHARACTERISTICS  METHOD  AND  ITS  APPLICATION 
TO  PRACTICAL  NONLINEAR  SUPERSONIC  PROBLEMS 
By  Antonio  Ferri 


SUMMARY 


The  method  of  characteristics  has  heen  linearized  hy  assuming  that 
the  flov  field  can  he  represented  as  a  basic  flov  field  determined  hy 
nonlinearized  methods  and  a  linearized  superposed  flow  field  that 
accounts  for  small  changes  of  boundary  conditions.  The  method  has  heen 
applied  to  two-dimensional  rotational  flow  where  the  basic  flow  is 
potential  flow  and  to  axially  synnnetrlc  problems  where  conical  flows 
have  been  used  as  the  basic  flows.  In  both  cases  the  method  allows  the 
determination  of  the  flow  field  to  be  simplified  and  the  numerical  work 
to  be  reduced  to  a  few  calculations.  The  calculation  of  axially  sym¬ 
metric  flow  can  be  siii5)lified  if  tabulated  values  of  some  coefficients 
of  the  conical  flow  are  obtained.  The  method  has  also  been  applied  to 
slender  bodies  without  symmetry  and  to  some  three-dimensional  wing 
problems  where  two-dimensional  flow  can  be  used  as  the  basic  flow.  Both 
problems  were  unsolved  before  in  the  approximation  of  nonlinearized  flow. 


INTRODUCTION 


The  use  of  the  method  of  characteristics  for  the  solution  of 
supersonic- flow  problems  requires  numerical  procedures  which  are  lengthy 
and  involved  and  which  must  be  repeated  for  each  set  of  boundary  condi¬ 
tions.  The  method  has  received  general  practical  application  only  for 
two-dimensional  or  axially  symmetrical  problems  in  steady  flow  and  one- 
dimensional  or  quasi- one-dimensional  nonsteady  flow^  and  only  very  few 
cases  of  general  three-dimensional  flow  have  yet  been  investigated. 

Many  problems  have  been  investigated  at  present  by  means  of  the 
linearized  theory  in  which  the  disturbance-velocity  components  (defined 
as  the  difference  between  the  local  and  the  free-stream  components  of 
the  velocity)  are  considered  and  are  small,  so  that  terms  of  second 
order  or  higher  can  be  neglected.  In  the  present  paper  a  simplification 
is  introduced  in  the  equations  of  motion  based  on  the  assumption  that 
one  of  the  velocity  components  or  the  variation  of  the  velocity  com¬ 
ponents  as  a  function  of  a  given  parameter  can  be  considered  small,  so 
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that  terms  of  second  or  hi^er  order  in  the  quantities  considered  small 
ceuti  he  neglected.  When  one  of  the  velocity  components -is  assumed  to  he 
small^  the  other  two,  velocity  components  can  he  esgpressed  in  two  parts, 
one  of  which  is  large  and  is  a  function  only  of  two  coordinate  positions, 
and  the  other  of  which  is  small,  of  the  same  order  as  the  third  velocity 
conqjonent,  and  is  a  function  of  all  three  coordinates. 

If  the  variations  of  velocity  coiaponents  as  functions  of  a  given 
parameter  are  considered  small,  all  three  velocity  con5)onents  can  he 
e3Q)re88ed  in  two  parts.  One,  large,  is  independent  of  the  parameter 
considered,  euid  the  other,  small,  is  a  function  of  the  parameter  con¬ 
sidered.  When  the  velocity  components  are  substituted  into  the  differ¬ 
ential  equations,  the  equations  can  he  divided  into  two  parts,  the 
differential  equations  'containing  the  velocity  conqjonents  considered  small 
become  linear;  therefore,  superposition  of  solutions  is  possible.  With 
this  assuiiq)tion  the  flow  field  can  he  represented  for  any  condition  by 
the  superposition  on  a  nonlinear  basic  flow  field  of  a  linearized  flow 
perturbation.  The  flow. field,  which  represents  the  variation  of  the 
basic  flow  due  to  the  changes  of  the  geometrical  parameter  considered, 
changes  linearly  with  the  parameter.  Because  of  the  simplification, 
the  superposed  flow  field  is  defined  by  differential  equations  of  hyper¬ 
bolic  type  Tdilch  have  characteristic  surfaces  equal  to  the  characteristic 
surfaces  of  the  basic  flow  field  and  known  coefficient;  therefore,  the 
superposed  flow  field  can  be  obtained  directly  without  the  iteration 
process  along  the  characteristic  net  of  the  basic  flow. 

A  particular  application  of  the  linearized  characteristics  method 
has  been  discussed  in  references  1,  2,  eind  3  in  which  bodies  of  revo¬ 
lution  at  small  angles  of  attack  have  been  considered.  In  the  present 
paper  the  basic  concept  of  the  linearization  is  discussed  and  exanqjles 
of  application  to  two-dimensional  rotational  flow,  to  conical  flow,  to 
axially  symmetric  flow,  and  to  some  general  three-dimensional  problems, 
are  discussed.  From  these  exanqjles,  other  applications  of  the  same  method 
to  supersonic  steady-  or  nonsteady-flow  problems  can.  be  visualized.  For 
exaii5)le,  the  method  can  also  be  applied  to  the  determination  of  the  flow 
field  in  supersonic  coiiq)ressors  or  turbines  having  supersonic  relative 
velocity  inside  the  passage.  In  this  case,  the  two-dimensional  flow  of 
the  cascade  of  the  blades  at  each  radial  station  or  the  axially  sym¬ 
metric  flow  can  be  assumed  as  the  basic  flow.  In  the  first  case  the 
radial  contponent  of  the  velocity  must  be  assumed  to  be  small,  while  in 
the  second  case,  the  tangential  component  of  the  velocity  must  be 
assumed  to  be  small. 


SYMBOLS 


a,b,c,d, e  coefficients  of  the  variables  in  the  characteristic  equa¬ 
tions  (defined  case  by  case) 
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n 


'll 


a  speed  of  sound 

An  coefficient  (  ^0%  +  ^o'^xSj 

B  and  C  coefficients  defined  ty  equations  (20)  and  (25) 

Cp  pressure  coefficient 

M  Mach  number 

R  gas  constant  or  radius  of  hodograph  diagram 

S  entropy 

u^v^v  velocity  con^onents  along  the  x-^  y-^  and  z-axes  in  Cartesian 

coordinates  or  along  x-  and  y-axes  and  peip)endicular  to 
the  meridian  xy-plane  in  cylindrical  coordinates 

Vr^Vn^'w  velocity  con^onents  in  polar  coordinates 

Vrp^v^  tangential  and  normal  velocity  conponents  in  front  of  the 

^  ^  shock 

V  intensity  of  the  velocity  vector 

\ 

P  Mach  angle 

7  ratio  of  specific  heats 

2  inclination  of  the  characteristic  lines  in  the  plane 

^  z  =  Constant  or  0  =  Constant 

G  coordinate  of  the  meridian  xy-plane  in  cylindrical  coordi¬ 
nates  or  of  the  rvl/’-plane  in  polar  coordinates 

*^^0  polar  coordinates 

cp  inclination  of  the  velocity  vector  with  respect  to  the 

x-axis 

components  of  the  rotation  along  the  x-_,  y-,  and  z-axes 

Subscripts: 

o  properties  of  the  basic  flov 

l,,.n  properties  of  the  superposed  linearized  flov  fields 


MCA  TN  2515 


TKE  EQUATIONS  OF  THE  LIHEAKI2ED  CHARACTERISTIC  SYSTEM 


Consider,  for  example,  a  flow  field  defined  by  a  velocity  vector 
V(u,v,w)  tbe  con^jonents  of  vhich  can  be  expressed  in  the  form 


J.1 

=  >"0 


V  =  v„  + 


o  +  /  ^n^n 
1 


^  ~  ^  ®n% 

1 

The  relation  between  entropy  and  rotation  states 


curl  V  X  V  =  grad  S 
yR 


Assiime  that  and  v^  are  functions  only  of  x  and  y,  then  by 

neglecting  terms  of  the  order  of  eqviation  (2)  becomes 


a^  _  „  /Svo  Suo 


Bx  yR 


PX  hy 


Svn  Sun\  (hvQ  Suo 


Sx  hy  )  \Sx  'by 


anvn 


M  a£  =  u  ( ^  ^ 

by  yR  °\bx  by 


M  a!  ^  ^  y-  .n  S~  a 

yR  °  2. —  \^y  bz  )  °  Z Sx  ) 
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Therefore^  the  entropy  field  can  he  expressed  as- 


Sx  Sx  I- —  ^  Sx 


Sy 


Sy 


>  (M 


Assume  that  each  coefficient  a^  .  .  .  is  constant  in  the 
entire  region  of  the  flow  field  where  it  is  not  zero  and  is  srnall^  so 
that  terms  of  the  order  of  a^^  or  higher  can  he  neglected.  In  this 

flow  a  basic  flow  field  exists,  represented  hy  the  velocity 
vector  Vq^Uq,Vq^  and  hy  the  entropy  distribution  on  which  a 

linearized  flow  is  supeiposed,  represented  hy  a  summation  of  N  three- 
dimensional  flow  fields,  each  of  which  is  proportional  to  the  corre¬ 
sponding  coefficient  a^.  The  basic  flow  field  is  a  two-dimensional 

flow  if  Cartesian  coordinates  are  used  or  an  axially  symmetric  flow  if 
cylindrical  coordinates  are.  used  and  can  he  determined  hy  known  methods. 
In  a  similar  way,  a  general  three-dimensional  flow  field  can  he  assumed 
for  basic  flow  if  the  flow  field  can  he  obtained  hy  sinple  analysis. 


The  equation  of  motion  obtained  from  continuity,  momentum,  and 
energy  equations  can  he  expressed  in  Cartesian  coordinates  in  the  form 


^(1  _  +  ^(1 . 

Sx  \  a2  y  Sy  \  a' 


5z 


a' 


uvlhu.  ^ 
a^ySy 


UWjSu  ^  ^w\  _  wv/5v  ^  _  Q 

a2\3z  hx)  a^ySz  Sy/ 


(5) 
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vhile^  in  cylindrical  coordinates,  the  equation  hecomes 


Sx 


w£  \  uv/5u  . 
a^y  bx) 


W_^  +  +  ^)  +  V  ^  0 

a^\y  dcp  dx  a2w  acp  dy  y 


(6) 


By  use  of  equations  (l),  expression  (5)  becomes 


where 


n 


_  7  -  1 
„  2 


+  ^o^n) 


The  basic  flow  is  a  two-dimensional  flow,  and  the  first  part  of 
equation  (7)  must  he  zero;  therefore,  the  second  part  of  the  equation 
must  also  he  equal  to  zero  and,  for  each  value  of  n. 
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^  4.  ^l^  -  +  — 

Sz  hx.  \  a.Q^  j  by 


+  .  j:o 


ao2/  ao2  Ux  Sy 


Suo/2uoUn  +  Uo?An\  ^  Svq /2voVn  +  Vp^An 


\  \  2 


by  bx  / 


Equation  (8)  is  a  differential  equation  of  hyperbolic  type  because 
Uq^  +  and  the  characteristic  surfaces  are  cylindrical  surfaces 

perpendicular  to  the  plane  z  =  Constant  vith  generatrices  coincident 
vith  the  characteristic  lines  of  the  two-dimensional  flow.  This  can  be 

seen  from  the  fact  that  the  coefficients  of  the  derivatives 

ox  dy 

Sv^  Su^ 

and  — ^  of  equation  (o)  are  expressed  as  functions  of  the  prop- 

Sx  8y 

erties  of  the  basic  flow  and  are  the  same  as  the  coefficients  of  the 

derivatives  which  define  the  basic  flow, 

ox  oy  oy  ox 

hvy.  ,  . 

and  the  coefficient  of  is  one  (see,  for  example^  reference  4, 

oz 

page  282).  Therefore,  the  disturbance  flow  field  Vj^^u^,Vj^,Wj^)  can  be 

obtained  by  the  method  of  characteristics  by  moving  along  the  charac¬ 
teristics  surfaces  vhich  are  cylindrical.  The  characteristic  net,  -which 
for  the  general  case  must  be  drawn  in  spatial  coordinates,  can  be  drawn 
in  this  case  only  once  for  any  value  of  a  and  n  and  is  equal  to  the 
net  of  the  two-dimensional  flow.  Equation  (8)  can  be  transformed  for 

practical  use.  In  the  equation,  the  terms  in  UqVq,  ^^9  ^^9 

ox  dy  ox 


and  ^ 
oy 

field. 


are  known  terms  and  are  given  from  the  two-dimensional  flow 


Assume,  in  the  plane  xy,  polar  conponents  for  V  and  cp  for  the 
velocity,  defined  by 


Uq  =  Vq  cos  cp, 


=  Vq  sin  cp, 
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and,  because  hl^er-order  terms  are  neglected, 

%  =  Vn  cos  cpo  -  cp^Vg  sin  9^ 
Vn  =  Vn  sin  cp^  +  cos  cp^ 


(10) 


iK 

t 


vhere 


V  =  Vo  +  a^Vn 
<p  =  cpo  + 


1 


(11) 


If  n  is  the  normal  to  the  projection  of  the  streamline  in  the 

plane  z  =  Constant,  the  derivative  ^  is  expressed  hy 

dn 


Ss  _  5S  V  I  dS  u 

Sn  hx  V  Sy  V 


(12) 


By  use  of  equations  (4)  and  (l),  equation  (12)  becomes 


^  =  .fe 

hn  VSx 


^  aSnUo  H^nVn  ,  f^Sp  ,  V"  ^ 

^  ^  +  Ha^V^  W  ^  j 


^0  + 


(13) 


By  considering  only  the  lowest-order  terms,  equation  (13)  becomes 


hS  ^  r—  ^ 

Six  Six^  I _  ^n  5ru  2 _ 


-  ^(Vo  -  ^oVn)  +  I^Mo  -  '"oVn)^ 


(U) 
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- 


’^o^n 


2/dSo  uo  ,  SSo  Vo  \  _ 


o  \Bx  Vo  By  V, 


because  the  ri^t  term  contains  as  a  factor  the  variation  of  entropy 
along  the  streamline  of  the  basic  flow;  therefore. 


Bvp  Bu^  _  ap^ 

Bx  By  7EV  BnQ  I 


+  A 

^'/RVo  Bno 


where  is  the  normal  to  the  streamline  of  the  basic  flow  in  the 

plane  xy. 

After  several  transformations  the  following  equations  can  be 
obtained  (see  reference  2):  In  the  plane  z  =  Constant  along  the  char¬ 
acteristic  line  defined  by 


-  =  X2  =  tan(p^  +  cp^) 


the  following  equation  is  valid: 


<ivn  _1  sin  pp  tan  Pp  ^ 
dz  Vq  cos(Po  +  cpo)  Vq  dx 


dcp,^  sin^p  I  dS„ 

tan  p-  +  - j - !-e — ^  -1-  _n 

^  dx  cos^cpQ  +  Poj  7R  dn 


^nBl  +  r  ° 

*  o 


and  along  the  characteristic  line  of  the  second  family  defined  by 


^  =  ^2  =  tan((Po  -  Po) 
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the  following  equation  is  valid: 


dw^  1  sin  tan  1  ,  1  ,  dcp^  sinBp^  1  ^Sn 

dz  Yq  cos(9o  "  Po)  ^o  ^  cos(cPo  -  p^)  yE  dn^ 

'Pn®2+r^2  =  0  (19) 

*  o 


where  and  Cg  are  coefficient  functions  of  x,  y,  Y^, 

2^,  and  along  the  characteristic  lines  of  the  first 

o  \dx  J^\dx  ! 

and  second  families  at  each  point  and  are  independent  of 
and  Wj^,  and,  therefore,  can  he  determined  once  for  the  basic  flow  and 
used  for  any  kind  of  disturbance  flow  in  the  limits  of  the  approximation 
accepted.  The  coefficients  Bj^,  32,  Cq^,  and  Cg  are  given  by  the 
following  expressions: 


1 

sinVo  <130 

’  1 

cos(cpo  -  Po) 

cos  Po  cos(9o  +  Po) 

7R  dn^ 

_ 

VoW  4 

sin  Po 

(20a) 


B2  = 


1 _ 

cos  cos(cp^  -  Pq) 


sinVo  dS^ 
7E  dn^ 


Pq) 

VqW  sin  p^ 


(20b) 


'  ? 
tan^pQ  + 


7  -  1 


2  sin^p^  cos^Pq 


2  <^0  sin3p^  ^ 

cos(cPo  -  Po) 

/l  ,  7  -  1  \ 

7R  dno  cos(cp^  +  p^)  Vo 

Idx  j 

2  cos(q3o  +  Po)cos2Pq 

i  2  sin2pj 

(20c) 
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Co  =  ^a.n%  +  + 


^  Vo\dx  l^\  2  sin^pQ  cos^Pq 


_2  ^  1  (^^o\  cos(cpo  +  Pq)  1 

7R  diiQ  cos(q3Q  -  Pq)  Voldx  cos((f>Q  -  p^)  cos 


{%  -  Po)  cos^p^ 


1  + 


7-1' 

2  sIh^Pq 


(20d) 


where 


\dx  /i 


is  the  derivative  along  the  characteristic  of  the  first 


family  and  y  along  the  characteristic  of  the  second  family. 


In  order  to  determine  the  value  of  w^^  at  each  point  of  the  char¬ 
acteristic  net,  the  following  relations  can  be  used  in  the  plane 
z  =  Constant: 


fc  -  5w  §v 
Sy  Sz 


^  -  Su  Sw 


> 


(21) 


5v 

Sx 


^u 

Sy 


and,  if  s^  is  a  streamline  projection  in  the  plane  z  =  Constant, 


^w  u  ^  5w  V 

hx  Y  5y  V 


u  ,  V  .  av  av  as 

V  V  az  az  az  7Rv 


(22) 
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In  the  approximation  accepted^ 


_  Y~  „ 
Bz  Z 


as  _  V~  . 

dz  z _  ^  az 


therefore^ 

i_  ^ 
'^o  S=o 


■where  is  the  streamline  of  the  hasic  flow. 

Equation  (24)  permits  the  determination  of  the  value  of  w  at  each 
point  of  the  plane  z  =  Constant  as  a  function  of  the  local  variation 
of  in  the  z-d±rection.  Therefore^  by  means  of  a  step-hy-step  pro¬ 

cedure^  all  the  flow  field  can  he  determined  hy  working  in  planes 
z  =  Constant  where  only  one  characteristic  net  must  he  used^  hy  means 

of  equations  (17)  and  (I9)  and  hy  determining  the  value  of  at  each 

point  of  the  net  hy  means  of  equation  (24).  The  calculations  are  started 
from  the  flow  field  defined^  or  along  a  surface  -which  is  not  a  charac¬ 
teristic  surface^  or  along  a  characteristic  surface  and  a  stream  surface. 
The  flow  field  at  the  starting  surface  must  he  determined  from  the 
boundary  conditions.  If  the  starting  surface  is  a  shock  wave,  the  flow 
at  the  shock  surface  must,  he  obtained  from  the  physical  properties  of 
the  shock  wave  related  to  the  boundary  conditions  considered.  Relations 
between  boundary  conditions  and  shock  waves  are  presented  in  detail 
subsequently  for  the  problems  considered. 

Similar  equations  can  be  obtained  by  using  cylindrical  coordinates 
in  place  of  Cartesian  coordinates. 

Equations  (l7)  and  (I9)  remain  the  same.  Only  the  first  term 

changes,  ^  becomes  and  the  coefficients  Bn  and  Bo  become 

oz  y  00  ±  d 


1  ^^n  _j_  ^^n  ^o 

Vq  Sz  Sz 


(24) 
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B 


I  _ 


sin  Po  sin(cpQ  -  p^) 


^  cos  Pq  cos(cpo  +  Pq)  y 


Bo‘  = 


sin  pQ  sln(cpQ  +  p^) 


+  B. 


COB  P^  cos(cp^  -  p^)  y 


(25) 


Equation  (25)  for  cylindrical  coordinates  becomes  (see  reference  2^ 
equation  (9c)) 


y  he  /Wq  y  Be 


Vp  sin  e^ 

y 


(26) 


The  use  of  the  linearized  characteristic  system  is  sinq)ler  in  many 
cases  than  the  complete  characteristic  system  and  reduces  in  some  prac¬ 
tical  applications  the  extent  of  the  numerical  vork  required,  especially 
if  solutions  of  many  similar  problems  are  required.  The  same  concept 
can  also  be  easily  applied  in  the  field  of  nonsteady  or  relative  motion 
of  flovs. 

The  flow  field  around  slender  bodies  without  axial  symmetry  can  be 
obtained  by  mjans  of  linearized  methods,  and  some  practical  three- 
dimensional  problems  not  analyzed  before  in  the  approximation  of  non- 
llnearized  flow  can  be  analyzed  by  this  method. 

In  the  next  sections  some  typical  possible  applications  are 
presented. 


SCME  TYPICAL  APPLICATIONS  OF  THE  LINEARIZED  CHARACTERISTICS  METHOD 
Two-Dimensional  Rotational  Flow  Fields 


Two-dimensional  potential  flow  permits  hodograph  solutions,  and, 
therefore,  any  kind  of  two-dimensional  supersonic  potential- flow  solu¬ 
tion  can  be  obtained  in  the  hodograph  plane;  the  numerical  solution  in 
the  physical  plane  then  requires  only  the  construction  of  a  character¬ 
istic  net  in  order  to  find  the  position  in  the  physical  plane  of  each 
^  poiJit  of  the  hodograph  plane.  The  solutions  of  problems  in  which 
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boundary  conditions  are  given  only  along  a  streamline  are  very  simple 

in  the  approximation  of  potential  flow  because  in  this  case  the  velocity 

is  constant  along  a' family  of  characteristic  lines  (single- wave  flow);  t 

however,  similar  calculations  for  rotational  flow  are  much  more  involved, 

because  a  step-by-step  procedure  is  required  for  the  solution  on  the 

holograph  plane  as  well  as  for  the  construction  of  the  characteristic 

net.  By  means  of  the  linearized  characteristics  method  a  rotational 

flow  field  can  be  considered  as  a  modification  of  a  potential  flow 

field,  and  the  linearized  superposed  flow  is  the  flow  which  takes  into 

account  the  effect  of  the  presence  of  shock  waves  and  the  effect  of 

rotation  in  the  flow. 

Consider,  for  example,  a  two-dimensional  profile  which  produces  a 
shock  wave  at  the  leading  edge  (fig.  l).  If  the  profile  is  curved,  the 
shock  is  curved,  and  the  flow  behind  the  shock  is  rotational.  Assume 
that  the  flow  field  behind  the  shock  can  be  expressed  as 


or 


^  =  '^o  +/  ^1 


V  =  v^  +  V- 


and 


^  =  ^0+^1  ^  9  =  q>o  +  9i  >  (27) 


where  u^  and  Vq  are  the  velocity  conponents  defined  by  a  potential 
flow  field,  Tidiich  in  this  case  is  a  single-wave  flow,  along  the 

characteristic  lines  of  the  first  family  is  zero,  and  and 

represent  the  flow  field  that  takes  into  account  the  reflections 
occurring  at  the  shock  and  the  effects  of  the  entropy  gradient.  The 
flow  field  represented  by  and  cp^  can  be  immediately  determined, 

and  Vq  and  cp^  are  constant  along  the  characteristic  lines  of  the 

first  family,  which  are  strai^t  lines.  The  characteristic  net  of  the 
flow  VqCPq  can  be  drawn  in  a  short  time.  Then,  along  the  characteristic 
of  the  first  family. 


sin^P 


tan  Pq  dcpi  +  ■  °  dS^  +  dx^^  +  ^  dx^^  =  0  (28) 
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where 


and 


Bl  =  - 


cos  pQ  cos(cpQ  +  Pq) 


U^\  -  Pq) 

yo\^  A 


2  sin  Pq 


=  _L  cos((Po  -  Pq)  _i _ A  ^  7  -  1  \ 

Vdx  /g  Vo  coB(cpo  +  Po)  cofl2p^^  2  Bin%l 


(29) 


(30) 


and,  along  the  characteristic  of  the  second  family, 

i  dV^  +  tan 
Vq  f 


sin^P^  Vi  dV„/  y  -  1  \ 

Po  ^'Pl  +  ‘^l  +  Po  +  - 1 - 1  “  ° 

''o  ’0  ^  2  sin^Pg  cos^P^y 


(31) 


Equations  (29)  and  (3I)  can  he  slii5)llfied  hy  introducing  the  gra¬ 
dients  (Vo)jj  and  ^cpQ  +  Pq)^  along  the  x-axis  defined  as 


-  (^L 


(■Po  +  Po)v  ■ 


y=0 

^(q^o  +  Pq) 

hx 


(32) 


y=0 


It  can  he  seen  from  figure  1  that,  at  any  point  A  of  ordinate  y^  on 
the  characteristic  line  a,  crossing  the  axis  at  A^,  the  vari- 

a 


at  ion 
is 


along  the  characteristic  line  of  the  second  family  at  A 


dx 


2Ja 


“  (Vox), 


1  - 


cot(cpo  +  Po) 

^o 

1 

cot(cPo  -  Po) 

^A 

1  _  _ 

(33) 


sin2(cpo  +  p^)^^ 
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i6 


MCA  TN  2515 


because  3q  and  qpQ  are  the  same  at  A  and  at  Aq.  Therefore, 


a 

1  - 


t 

(34) 


■vdiere  a  and  b  are  constants  along  each  characteristic  line  of  the 
first  family.  Then  equation  (28)  along  the  characteristic  line  of  the 
first  family  becomes 


tan  Pq  dcp^  + 


sln^pQ 

7R 


dS2_  -  cp-]^ 


cos  (cpQ  -  Po)  _ _  a  dy  ^ 

sin(cpo  +  Po)  pQ  Pq  ^  “  by 


hfl  +  y  -  1  \  a  dy  cos(cpo  -  pp)  q, 
^o\  2  sin^Po/b  -  by  sin(cpQ  +  Pq)  cos^Pq 

Therefore^  if 


a 


and 


_  foo 

sin2(cpQ  +  p^) 


(35) 


(36) 


the  equation  along  the  first  characteristic  line  becomes 


^  -  tan  pQ  dcpi  +  dS  +  cp  ^  d  log(l  -  by) 

Vo  7K  ^  X  D 


1  +  - }  \ta.n  pQ  ^  -  d  log(l  -  by)  =  0 

2  sin  p„/  ''o 


(37) 


>• 


S’- 


3 


NACA  TN  2515 


17 


and^  along  the  second  characteristic  line, 


— ^  +  tan  P  dcp-  + 
V  o  1 

^o 


sin^p^  Vn  dV^  /  7-1 

- ^  dSn  +  — ^/tan^P^  +  — - ^ - r — 

^o  ^o  I  2  sin^pQ  cos^Pq 


All  the  coefficients  of  equations  (37)  3*nd  (38)  constant  along 

characteristic  lines  of  the  first  family  and  can  he  calculated  at  few 
points  on  the  x-axis. 

The  coefficients  of  equations  (37)  (38)  independent  of  cp^ 

and  therefore,  with  one  calculation  from  a  point  A  and  C  of  the 

net  of  figure  1  the  values  of  and  cp^  at  a  point  B  can  he 

obtained  directly  without  the  necessity  of  an  iteration  process, 

p 

including  also  terms  of  the  order  of  (Ax)  in  each  step.  Indeed,  if 
all  the  quantities  which  are  variable  along  the  characteristic  lines  are 
expressed  in  the  form 


"3  ■  “a  +  (Ij).  ^  S  ^  *  °(**^) 


(HI  ■  ('* 


Ax  +  o(ax^) 

/a.  Sx^ 


then 


“B  ■  ^A  '  A*  *  to  I  +  o(ix3)  (40 

Consider  now  equations  (37)  and  (38).  At  point  C  equation  (37)  has  the 
form 


^^1)  a  +  1^^  h  +  0)  +  V 

8x  /c  \8x  /q  C  IcW 


4 
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•while  at  point  B  it  has  the  form 


then,  from  equation  (40), 


Indeed, 


r 


* 
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(1^2) 

and  equation  (38)  along  the  second  characteristic  line^  terms  of  the 
order  of  (xg  -  being  neglected,  can  he  written  in  the  form 


20 
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Because  of  the  possibility  of  considering  directly  in  the  calculations 

terms  of  the  order  of  (Ax)^  also,  large  steps  can  be  used  in  the 
characteristic  net  and  the  effect  of  entropy  gradients  can  be  easily 
evaluated,  after  the  basic  characteristic  net  and  a  fev  streamlines  of 
the  basic  flow  have  been  determined  and  the  coefficients  of  equa¬ 
tions  (39)  and  (i)-0)  calculated  at  a  few  points  on  the  axis. 

All  the  coefficients  are  constant  along  each  characteristic  line 
of  the  first  family;  therefore,  in  going  from  B  to  D,  only  the  terms 
containing  and  y  must  be  changed  and  the  calculations  are  slinpli- 

fied  to  some  extent  with  respect  to  the  rotational- flow  characteristic 
calculations. 


Axially  Symmetric  Plow  Fields 

For  axially  symmetric  flow,  the  equation  of  linearized  character¬ 
istics  becomes; 


1  ^^1 


-:_2 


dcpi  sin^p  dSi  ^  ,  Vn 
~  r-i  +  cPjlEj^  + 


^ - —  -  tan  3^ - ,  - 

Vo  dx  dx  /R  dx  "1"1  ■  Vn  "1 


C  =  0 


m 


along  the  characteristic  line  Xj_  =  tan^P^  "Po) 


(1^5) 


along  the  line  Xg  =  ■tan(cpo  -  p^)  vhere  B^’,.  ,  C^,  and  Cg  are 

defined  by  equations  (25)  and  (20). 

The  introduction  of  equations  (44)  and  (45)  siiiq)lifies  noticeably 
the  numerical  calculations  without  affecting  sensibly  the  precision  of 
the  results.  The  practical  use  can  be  as  follows:  A  basic  body  shape 
la  determined  first  by  means  of  characteristic  calculations,  nnii  the 
characteristic  net  is  then  obtained.  The  basic  calculations  must  be 
extended  in  a  region  in  front  of  the  shock  wave  determined  with  the 
usual  procedure  (for  example,  reference  3),  and  in  a  region  inside  the 
body  as  shown  in  figure  2  in  order  to  determine  all  the  flow  field 
necessary.  Thus,  a  reduced  number  of  points  of  the  characteristic  net 
are  chosen  at  idilch  the  superposed  flow  field  for  each  boundary  condi¬ 
tion  different  from  the  basic  shape  will  be  determined.  The  number  of 
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points  required  depends  on  the  magnitude  of  the  supeiposed  flov  field; 
however,  the  number  is  usually  small,  because  for  each  step  between  two 
%  points  A  and  B  the  disturbance  velocity  can  be  expressed  in  the  form 


ential  equations  (kk)  and  (^5)^  the  coefficients  of  the  differential 
equations  are  independent  of  the  solution  and  are  known  at  both  points. 
If  the  entropy  term  is  neglected,  by  applying  finite-difference 

methods  the  velocity  components  at  a  given  point  g5  of  the  character¬ 
istic  net  can  be  expressed  from  the  values  at  two  points  f5  and  gh  in 
the  form: 


(''1),,-  (''i 

ary  conditions  considered,  "while  the  coefficients  L,  m,  n,  p,  q, 
and  r  are  functions  only  of  the  basic  flow  field  and,  therefore,  must 
be  determined  only  once  for  any  boundary  condition  considered.  These 
coefficients  are: 


where  (cp^) 


)  are  functions  of  the  hound- 
^f5 
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aD-|  ,  . 

If  the  entropy  tenns  in  are  considered,  two  more  terms  in  (S;|_)^^ 

and  m.Bt  he  considered  in  equations  (47)  and  (48)  and  two  more 

terms  in  equation  (49),  However,  because  these 

terms  are  small,  they  can  usually  be  neglected  in  practical  calculations. 


Because  all  the  coefficients  are  determined  only  once,  the  deter¬ 
mination  of  or  (^1)^1^  each  boundaiy  condition  is  simple. 

Points  at  the  boundary  can  be  investigated  by  means  of  equations  (47), 
(48),  and  (49).  Points  on  the  shock  can  be  analyzed  in  a  similar 
manner.  A  practical  calculation  can  be  performed  in  the  following  way: 
The  basic  flow  field  and  the  number  of  points  in  which  the  superposed 
flow  field  will  be  considered  having  been  determined,  the  new  boundary 
conditions  (shape  of  the  body)  are  placed  in  the  characteristic  net 
(fig.  3)*  If  A  is  the  point  where  the  basic  body  departs  from  conical 
shape  and  OBC  is  the  new  boundary  condition,  the  flow  field  between  the 
surface  of  the  conical  body  OB  and  the  conical  shock  OE  are  known  from 
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cone  calculations;  therefore,  the  values  of  V  and  cp  at  each  point  Ba, 
la,  2a,  'and  Ea  of  the  basic  characteristic  net  are  determined,  and 
by  difference  the  values  of  and  can  be  obtained.  From  al  the 

flov  at  the  point  L  of  the  boundaiy  can  be  determined  from  equation  (48) 
where  cp^  is  known  at  L,  and  the  coefficients  at  L  can  be  determined  by 

linear  interpolation  between  bl  and  cl.  From  the  values  of  and  cp-j^ 

at  al  and  LI,  the  corresponding  values  at  bl  are  interpolated.  Then  all 
the  values  for  the  line  b  can  be  obtained.  For  the  determination  of 

and  cp^  at  a  point  F  on  the  shock  wave,  the  equations  of  the  shock  wave 

and  equation  (47)  are  valid.  At  the  point  F  the  values  of  cp^  and  Vq 

are  known,  and  from  the  equations  of  the  shock  wave  the  value  of  —  as 

8cp 

function  of  cp  can  be  determined.  (The  value  of  cp  fixes  the  deviation 
across  the  shock  wave.)  Therefore, 


(V)r  -  (V„)^  +  (Vi)^  =  (v„)^  +  (|)(T0. 


Then,  equation  (47)^  applied  between  the  points  4b  and  F,  gives 


=  Hi  -  ("i)h  Hi  -  (%).  ft^ 


and  the  value  of  can  be  determined. 

The  work  required  in  the  calculation  of  the  flow  field  for  the 
basic  body  and  the  determination  of  the  coefficients  L,  m,  n,  p,  q, 
and  r  can  be  reduced  to  a  minimum  if  conical  bodies  are  assumed  as 
basic  bodies  for  the  calculations,  because  in  this  case  the  basic  flow 
is  available  in  tabulated  values  (reference  5)  'the  coefficients  L, 
m,  n,  p,  q,  and  r  are  functions  only  of  the  polar  coordinate  \lr. 

For  conical  flow,  the  calculations  can  be  performed  in  the 
following  way:  From  coni  cal- flow  calculations,  the  values  of  v^  and 

Vj^  as  functions  of  ^  are  known,  where  v^  is  the  radial  coii5)onent 
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9^*3.  "tlie  n.oz'ina.l  con^onen't  of  tlio  flow  field  in  polar  coordinates 

referred  to  the  limiting  velocity  (fig.  4).  The  following  expressions 
^o>  *^0  determined  from  the  cone  calculations: 


V2  = 


+  V, 


2  ~  3-  +  y  ^  sin  Pq 

Vo 

=  (p  +  7 


'  >  (53) 


tan  7  =  ^ 
^r 


Therefore,  +  9^  and  9^  -  are  known  as  functions  of  f.  From 

a  point  A  on  the  conical  body,  the  characteristic  line  of  the  first 
family  AC  and  of  the  second  family  AD  can  be  drawn.  The  line  is  defined 
by  the  expressions 


■  ^A  _  +  *p)e 

*E  ”  ^A  2 


Vj  -  -  'p)a  +  "  'pV 

2 


In  order  to  construct  a  characteristic  net  that  requires  only  a  small 
amount  of  calculation,  the  points  of  the  net  are  chosen  along  strai^t 
lines  from  0  so  that  the  conical  property  of  the  flow  can  be  utilized. 
The  net  can  be  constructed  by  fixing  the  steps  along  the  body. 

When  point  B  is  chosen  along  AO,  the  point  E  along  AC  is  determined 
by  drawing  BE  parallel  to  AF.  From  E  and  F  (along  OE  and  AD)  the  lines 
GE  parallel  to  LF  and  FG  parallel  to  EH  can  be  determined,  and  the 
point  G  can  be  obtained..  From  G  and  H,  point  M  can  be  obtained,  and  by 
proceeding  in  a  similar  way,  all  the  characteristic  net  can  be  deter¬ 
mined.  The  coefficients  of  equations  (50)  are  the  same  for  each  point  N, 
F,  and  E  along  the  same  radius.^  The  tabulated  values  are  given  only  in 
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the  region  between  the  body  and  the  shock  wave;  however,  the  calcu¬ 
lations  can  be  extended  by  means  of  conical  calculations.  For  example, 
i  from  the  following  eq.uations  (see,  for  exan^jle,  reference  k,  p.  243, 

and  the  following  pages) 


and 


(R) 


Vr  +  co't  i 


1  - 


2v. 


n 


(7  -  l)(l  -  Vo^) 


(54) 


tc 


y  (55) 


At  -  (E  -  CO,  At  .  (E)^^ 


(where  R  is  the  radius  of  curvature  of  the  streamline  in  the  holo¬ 
graph  pl6Lne);  therefore,  the  characteristic  net  can  he  extended  to  the 
outside  flow.  For  conical  flow  the  coefficients  B2%  and  C2 

can  he  determined  as  functions  of  For  conical  flow,  coef¬ 

ficient  Bq_*  of  equation  (25)  becomes 


Bl^ 


1 

sin  sin((Po  -  p^) 

cos  3^  cos(9„  +  p^) 

y 

i  ^°K^0  -  Pq) 

Vo  dt  Vdx/^^  sin 
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and  C2_ 


"but 


and 


Therefore 


sin  Ilf  sin^cpQ  +  -  i|f) 

y  cos(cPo  +  Po) 


sin  sin^’l"  -  cp^  + 

y  cos(cpo  -  Pq) 


(56) 


(57) 


1  ^  sin  Pq  sin(q)o  -  pp)  ^ 

cos  Pq  sin  \|f  sin^cpQ  +  Pq  - 


l'^ 

1  1 

-  CPo  +  Pol 

\dt  j 

i|f  sin  Pq  sin 

(cpo  +  Pq  -  Ilf) 

(58a) 
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and 


Vj.  +  cot  ilf 


(7  +  l)(l  -  Vo^) 


4 


If  equations  (58)  are  used,  the  coefficients  of  equations  (50)  cap  be 
e:^ressed  as  follows: 


"  -  -  M  (v  *  “gs) 

‘^g5) 

"  J('''g5  ‘  '''f5)  {^f5  ^g5) 

""  "  i  (''g5  "  ''f5)  ('^f5  ^5) 


>  (59) 


Because  the  coefficients  a,  b,  c,  and  d  are  functions  only 
of  lie  and  of  the  free- stream  Mach  number  M,  they  need  to  be  calculated 
only  once  for  different  values  of  and  given  in  tabulated  form,  and, 

therefore,  the  calculation  of  any  flow  field  for  vdiich  the  basic  flow 
can  be  considered  a  conical  flow  can  be  reduced  to  the  solution  of  a 
few  linear  equations  with  known  coefficients. 


In  order  to  try  the  method,  the  flow  around  an  ogive,  as  shown  in 
figure  5^  has  been  determined  at  M  =  3 -016,  by  the  method  of  character¬ 
istics  and  by  the  method  of  linearized  characteristics.  In  figure  5  the 
usual  characteristic  net  is  shown,  while  in  figure  6  the  linearized 
characteristic  net  and  the  basic  body  are  presented.  As  a  first  basic 
body  the  cone  chosen  is  the  cone  tangent  to  the  apex  having  a  cone  angle 
of  12.5°.  The  cone  chosen  is  not  the  most  convenient  because  the  values 
of  at  the  end  of  the  ogive  are  large,  and  a  better  approximation 

would  be  obtained  if  a  cone  of  smaller  cone  angle  would  be  considered  as 
the  basic  body.  A  cone  some'vdiat  different  from  the  ogive  considered  has 
been  chosen  in  order  to  have  some  information  on  the  approximation  of 
the  method  for  sensible  variations  of  the  shape  of  the  body  from  the 
basic  body. 
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Downstream  of  the  characteristic  line  AB,  the  velocity  com¬ 
ponents  Vj  and  cp^  would  become  large  because  at  the  surface  of  the 

body  the  component  cp  is  quite  different  from  the  con5>onent  cp^  of 

the  basic  flowj  therefore,  the  flow  determination  in  the  region  down¬ 
stream  of  the  line  AB  has  been  considered  as  a  new  problem,  defined  by 
the  flow  along  AB  and  from  the  streamline  that  represents  the  body 
shape.  In  this  region  a  new  basic  flow  has  been  considered.  Again,  a 
conical  flow  field  has  been  assumed  as  the  basic  flow.  From  the  values 
of  <p  and  V  at  B  and  A,  and  from  the  order  of  magnitude  of  cp  down¬ 
stream  along  the  body,  a  conical  flow  field  that  would  give  a  small 
disturbance  component  in  this  region  has  been  selected.  The  cone  chosen 
for  the  second  part  is  a  5°  cone  at  M  =  3.077.  The  cone  is  entirely 
contained  within  the  body  considered  in  the  region  used  in  the  calcu¬ 
lations  as  shown  in  figure  6,  and  the  most  convenient  region  of  the 
conical  flow  field  is  used  for  the  calculations. 

In  order  to  pass  from  one  basic  body  shape  to  the  other,  the  com¬ 
ponents  and.  cp^  along  the  characteristic  line  AB  for  the  second 

basic  body  must  be  determined.  This  operation  can  be  performed  by 
determining  the  value  of  Vq*  +  7’]_'  and  <Pq'  +  cp]^’  at  the  points  a, 
b,  and  n  of  figure  6  for  the  first  basic  body,  by  Interpolating  the 
values  at  the  points  a*,  b',  and  n’  between  characteristics  1  and  2 
along  the  characteristic  of  the  other  family  and  then  determining  the 
new  value  of  and  cp^"  for  the  second  basic  body  along  AB  from 

the  expressions 

=  V  -  Vq" 


9i"=cp-cPo" 

and 

V  =  Vo*  +  Vi* 

9  =  ^>0'  +  fPi' 

where  Vq"  and  cpo"  are  the  values  of  V  and  cp  for  the  second  basic 
body  at  the  position  considered  and  V^*  and  cp^*,  the  values  for  the 
first  iDasic  body. 
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The  calculations  by  linearized  characteristics  required  the  solu¬ 
tion  of  11  linear  equations  of  the  type  of  equation  (^7),  or  (i)-9) 

with  one  unknown  and  the  Interpolation  of  four  points,  which  can  he  done 
in  a  very  short  time  (of  the  order  of  1  hour)  when  the  net  is  drawn  and 
the  coefficients  of  the  equations  are  determined.  The  pressure  distri¬ 
bution  obtained  is  presented  in  figure  7  and  is  compared  with  the  pres¬ 
sure  distribution  obtained  from  the  exact  method. 

For  the  back  part  of  the  body  where  the  flow  differs  slightly  from 
parallel  flow,  a  cylindrical  body  with  uniform  flow  at  different  Mach 
numbers  can  be  considered  as  the  basic  body.  In  each  region  of  the 
flow  the  Mach  number  for  the  basic  flow  is  constant;  however,  the  com¬ 
ponent  u  can  be  approximated  conveniently  by  changing  the  basic-flow 
Mach  number. 

When  the  body  has  a  tail,  a  conical  solution  as  proposed  in  refer¬ 
ence  6  for  flow  Inside  a  tube  can  be  assumed  as  the  basic  body.  In  this 
case,  each  streamline  of  the  conical  solution  can  be  considered  as  the 
shape  of  the  basic  body  (fig.  8),  and  the  flow  field  can  be  obtained 
from  conical-flow  calculations,  which  can  be  determined  easily  in  the 
hodograph  plane.  By  changing  the  strength  of  the  final  shock  of  the 
conical  solution,  different  ratios  between  maximum  cross-sectional  area 
and  tail  area  can  be  obtained.  Any  part  of  the  streamline  can  be 
assumed  as  the  basic  body  shape. 

The  present  method  does  not  require  the  existence  of  a  linearized 
solution  and,  therefore,  can  be  applied  also  at  hl^  Mach  numbers. 

This  method  permits  obtaining  the  shape  of  the  shock  wave  and  taking 
into  account  entropy  variations.  High  precision  can  be  obtained  by 
using  several  basic  flow  fields  for  the  different  regions  of  the  body 
considered.  Because  of  the  sinqjlicity  of  the  calculations,  the  systemr- 
atic  calculations  and  tabulation  of  coefficients  of  equations  (58)  for 
different  cones  and  different  Mach  numbers  would  be  of  great  practical 
interest. 

Tabulated  values  can  be  obtained  also  in  the  following  way:  For 
each  cone  OC  of  cone  angle  cp^  considered,  a  superposed  flow  must  be 

calculated  as  shown,  for  example,  in  figure  9(a).  The  values  of 

and  cp^  for  this  superposed  flow  are  obtained  at  given  points  of  the 

characteristic  net.  Because  of  the  linearization  of  equations,  if  the 
superposed  flow  changes  in  Intensity,  the  values  of  and  cp^  at 

every  point  change  proportionately.  Because  of  the  conical  property, 
if  the  point  A  (fig.  9(b))  moves  along  OC,  the  flow  field  changes  in 
scale;  therefore,  the  effect  on  a  point  E  due  to  the  supeiposed  flow 
field  ZJfp^  starting  at  A  is  equal  to  the  effect  of  a  linearized  flow 
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field  starting  at  A*  and  of  Intensity  at  a  corresponding 

point  E'  defined  by 


V  =  V  ^  V 


Therefore,  when  the  flow  field  for  the  disturbance  AB  is  determined, 
the  effect  of  any  disturbance  of  the  type  of  the  disturbance  AB  in  the 
entire  flow  field  can  be  obtained.  Then  any  body  shape  can  be  con¬ 
sidered  as  a  superposition  of  flow  fields  of  the  type  of  flow  corres¬ 
ponding  to  the  disturbance  placed  along  the  basic  cone.  From  the 

simple  calculation  of  the  flow  for  the  shape  AB,  the  velocity  can  be 
determined  by  means  of  the  equation 


V  =  V, 


n 


n 


where  n  is  the  numher  of  the  superposed  flows  that  affect  the  point 
considered. 

^  Consider^  for  exaii5)le^  figure  10.  Several  linearized  flow  dis¬ 

turbances  must  he  superposed  on  the  basic  conical  flow  field.  First, 
a  superposed  conical  flow  at  0  that  can  be  obtained  from  conical- flow 
'  calculations  must  be  considered.  At  a  superposed  disturbance  must 

be  added  in  order  to  satisfy  the  boundary  conditions  at  A.  If  the 
calculations  have  been  performed  for  the  disturbance  at  a^  and 

for  =  1,  the  velocity  at  A  can  be  obtained  from 


’a  -  (Vi)^  t  (V2)^ 


Where  (Vq)^  is  the  velocity  of  the  basic  flow  at  A,  (V^)^  is  the 

velocity  disturbance  at  A  due  to  the  conical  flow  superposed  at  0, 
and  (^2)^  proportional  to 
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and  can  be  obtained  from 


iH 


vhere  b  is  determined  by 


In  a  similar  way  V2  is  determined  at  the  points  C,  and  so  forth. 
At  B  another  linearized  flow  field  having  velocity  coii5)onents 
and  cp^  must  be  considered  where 


(”3)^  ■ 


(’■0)3  +  "  ('^2)3 


and 


(’3)b  = 


(’’3)b 


where  b  is  defined  by 


"  yb 


Because  of  the  rapidity  of  calculation,  the  variation  of  any  geometrical 
parameter  can  be  investigated  in  practical  applications  without  the 
necessity  of  a  large  amount  of  numerical  work. 
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Conical  Flov  Field  without  Axial  Symmetry 

The  calculation  of  slender  bodies  without  axial  symmetry  requires 
the  determination  of  conical  flow  without  axial  symmetry^  which  can  be 
done  by  means  of  the  linearized  characteristic  method.  The  basic  prob¬ 
lem  of  the  determination  of  conical  flow  consists  in  determining  the 
shape  of  the  conical  shock  wave  produced  by  the  body.  When  the  shape 
of  the  shock  is  determined,  the  flow  field  around  the  body  can  be 
obtained  by  means  of  numerical  calculations  (see,  for  exaii5)le,  refer¬ 
ence  4).  .  Because  the  relation  between  the  shape  of  the  body  and  the 
shape  of  the  shock  wave  is  not  known  a  priori,  the  method  of  linearized 
characteristics  can  be  particularly  useful  for  flow  determination  of 
this  kind.  An  approximate  shape  of  the  shock  wave  is  assumed  as  the 
basic  solution  and  the  flow  field  inside  the  shock  is  determined;  then 
a  linearized  flow  field  is  supeipjosed  in  order  to  satisfy  the  boundary 
conditions  at  the  body.  The  calculations  are  siiip)le  if  the  basic  flow 
can  be  determined  analytically  or  numerically  without  a  large  amount  of 
calculation.  For  example,  for  slender  bodies  the  basic  flow  can  be  the 
axially  symmetric  flow  for  which  values  are  available  in  tabulated  form. 
Consider  a  conical  shock  wave  which  can  be  defined  in  polar  coordinates 
as 


00  00 

=  (^o)„  +  H  (^n)„  +  51  (^m) 


m0 


(60) 


where  all  the  values  of  and  are  small  so  that  terms  of  the 

order  of  can  be  neglected.  Such  a  shock  wave  is  approximately  of 

circular  cross  section,  as  is  found  for  slender  conical  bodies.  If  the 
flow  is  assumed  to  have  a  symmetry  plane,  the  second  summation  of  equa¬ 
tion  (60)  is  equal  to  zero. 

The  velocity  components  in  the  radial  direction  “the 

tangential  direction  direction  normal  to  the 

shock  front  of  the  shock  wave  are  (see  fig.  11):. 


cos 


sin  >lr  sin  a 


(61) 
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•where  all  the  velocity  components  are  referred  to  the  limiting  velocity 
and  a  is  the  angle  of  the  dihedral  het'ween  the  plane  nn-rmni  to  the 
shock  wave  and  the  plane  containing  the  reference  axis.  Across  the 
shock  wave  the  following  relations  are  valid; 

(''1)2 

^t)i 

(Ml 

If  the  velocity  con^onents  v^.,  v^^,  and  w  of  the  flov  field  in  polar 

coordinates  are  considered  at  each  point  ■'k  =  ’'Ic^  +  ^  nijr  cos  n6 
(fig.  11), 

=  ^r  =  Vi  cos  y  (53^) 

y  ^  _  7-1  (l“  Vj^^cos^  -  V^^sin^'f  sin^^ 

rn)^  -  -  TTT  - ^  sin  ilr - ^  ‘  ^ 
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Therefore^  the  velocity  components  behind  the  conical  shock  wave  of 
equation  (60),  if  terms  of  the  order  of  are  neglected,  are 
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Then,  the  velocity  components  hehind  the  shock  defined  by  equation  (60) 
can  he  expressed  as: 


O  T  / 


>  (65) 


(w)^^  =  ^  (t^)  sin  ne 


where  (v^)^,  (v^)^,  and  Wj^  are  independent  of  Q, 

v_  and  Vn  correspond  to  the  flow  field  for  circular  conical 
o  o 

shock  aITq  ,  and  ('v^)  ,  (vj^)  ,  and  (-w)  at  the  circular  cone  I'’  = 

n  '  'n  ^s 

are  independent  of  n  and  are  defined  by 
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It  can  be  shown  that  the  flow  field  defined  as  in  equation  (65)  satis¬ 
fies  the  boundary  at  a  surface  of  a  conical  body  defined  by 


in  =  ir.-  +  >  cos  ne 


if  the  terms  are  small  and  quantities  of  the  order  of  or 

hi^er  are  negligible  with  respect  to  terms  of  the  order  of  ilr  .  At 

n 

the  surface  of  the  body  the  velocity  must  be  tangent  to  the  surface  of 
the  body;  therefore,  at  each  point  A  of  the  body 


A  sin  d0 


sin  n0 
sin  tr. 


and,  therefore,  the  boundary  conditions  can  be  e3q)ressed  in  the  approxi- 
mation  considered  as 


W*  -  (X) 


t  \  at 
°c  °c 


r  (^n), 


cos  n0  + 


sin 

where  the  subscript  ilf^  indicates  quantities  at  the  surface  of  the 

basic  circular  body  and  indicates  quantities  at  the  surface  of  the 

conical  body  considered  where  the  parameters  (tn)^  are  given  by  equa¬ 
tion  (60)  and  define  the  shape  of  the  shock.  Because  each  term  of  the 
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ri^t-side  siinmation  is  of  the  order  of  boundary  conditions 

are  ' 


Equations  (68)  show  that  the  basic  flow  defined  by  the  components  v_ 

^o 

and  is  the  flow  corresponding  to  a  circular  cone  of  angle  iITq 

and  permits  relation  of  the  equation  of  the  conical  shock  to  the  equa¬ 
tion  of  the  conical  body.  The  coefficients  (v^.)  ,  (vj^)  ,  and  (w)^ 

for  different  values  of  n  can  be  determined  for  a  given  free- stream 
Mach  number  and  value  of  or  (tITq)  ;  therefore,  from  equa¬ 

tions  (60),  (65),  (67),  and  (68),  the  flow  field  around  any  conical  body 
of  the  type  given  by  equation  (67)  can  be 

are  small. 

The  determination  of  the  quantities 

functions  of  If  can  be  obtained  from  the 
ence  3) J 


obtained  when  the  terms 


following  equation  (see  refer- 


8v  / 

+  v_  cot  ilr  +  ^—2.11 
n  ^  hif\ 


Znf)  ^ 

a2  /  sin  ilr  d0  \ 


a2  ySi|r  sin  ^  h9 


(69) 
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Equation  (69),  because  of  expressions  (65),  gives 


H-n) 


a  i-:£a 


V^o  ao*^  /  Y  2 

-I  ^‘■o 

a  2 

i_  o  _J 


M  /  ,r-l^n„->Vn„cot»_  ; 

'  -  4  ..2  roX 


^O 


sin  ilf 


Because  of  the  approximation  considered,  the  entropy  remains  constant 
in  each  meridian  plane  outside  of  a  vorticose  layer  of  infinitesimal 
thickness  around  the  body  (see  reference  3).  Indeed, 


a£  ^ 
7R  Se 


^n  sin  t  -  v^  -  v^  ^  +  Vj.v  sin  ^  +  v^v  cos  t 


Vy.  sin  Hr  ^  =  -  v  ^ 

M  be 


From  equation  (71)  there  results  ||  of  the  order  of  ;  there- 

fore,  'where  Vj^  jf  0,  ^|j-  is  of  the  order  of  and  can  be  neglected. 

Then 


a  oo  _  A 

7R  ” 


TT  j.  w  ^n 


-  -w  —  + 


Si|r  sin  fhe  ^ 


+  v_v„  -  w2cot  ''|r 
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or,  in  the  approximation  considered, 


ov-n 

Then  Vj-  +  — | —  is  the  radius  of  the  hodograph  diagram  in  the  plane 
0  =  Constant  (see  reference  !)•)  and  and  (Vj.)  can  be  obtained 

from  a  step-by-step  calculation  from  t  to  '>1^0  =  ‘tn  by  means 

s  c 

of  the  equations 


vhere  at  station  ^  is  obtained  from  equation  (70)  and 


can  be 


calculated  from  the  values  of  (''^n)  >  ('”')n.  ^ 


The  value  of  ^  obtained  from 


— ^1  can  be  obtained  from  equation  (71)  vhere 

If 


and 
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S  -  Sq  +  S-L 


COS  n0 


and  is  independent  of  t  and  can  "be  detennined  from  the  equations 

of  the  shock  from  the  expression 


q  ^  /dS 

■  te; 


where  is  the  inclination  of  the  shock.  Then 


U  =  v_  sin  i  +  v~  (v-)  +  Vt,  (vti)  + 

1  ^o  ^o'  ^'n  “o'  “'n 


Vr^Wn  sin  >|r  +  Vn^w^  cos  i 


vhile^  if  the  quantity  is  neglected  and  the  flov  is  considered 

potential  flov^ 


(^r)n 


n  gj_n  ^ 


The  method  presented  has  been  applied  to  the  determination^  for  the  con¬ 
dition  of  zero  angle  of  attack^  of  the  flov  field  around  a  cone  having 
an  elliptical  cross  section  with  axes  in  the  ratio  of  about  1  to  3^  i*or 
which  experimental  data  were  available  at  M  =  1,8^  and  for  a  cone  of 
elliptical  cross  section  of  ratio  1  to  1.88.  The  calculations  have 
been  performed  in  the  following  way:  The  bodies  are  shown  in  fig¬ 
ures  12(a)  and  12(b).  The  value  of  ijr  at  0  =  0^  is  equal  to  6.30^ 
while  the  value  of  t  at  0  =  90^  is  equal  to  l8.4^. 

The  angle  ^  can  be  expressed  as 

=  \|r  +  'iIt-i  cos  20  +  io  cos  kO  +  cos  60  + 


tl,  cos  80  +  cos  100 


(81) 


k2 
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but 


tan  t  =  tan  + 


tan 


cos 


therefore,  by  determining  the  value  of  ijr  at  six  points,  the  following 

values  have  been  obtained:  If  i|fQ  =  10°, 

C 


=  -5.07° 

=  1.9i^° 

t,  =  -0.84° 

3c 

=  0.43° 

%  =  -0.16° 

5c 


Therefore,  a  10  cone  at  zero  angle  of  attack  is  assumed  as  the  basic 
body.  Eeference  5  gives  tabulated  values  for  a  10°  cone  at  =  I.816 

and  the  calculations  have  been  performed  at  this  Mach  number.  The 
table  gives 


%  =  34.45° 

°s 


The  entropy  variations  are  small  and  are  neglected.  Then, 


(w)  =  0.0255 

'  1-265 
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The  values  of  ('v„')  ,  (y„)  ,  and  (-w)^  between  \lf  and  have 

V  \  u  Og  Oq 

been  determined  by  means  of  equations  (TO),  (75),  and  (80).  Then  , 

s 

1^2^,  and  have  been  determined  from  equations  (68)  where 


-2  X  0.6000 


The  values  of  obtained  are: 


%  =  -0.24° 

-‘■s 

io  =  5.1°  X  10“3 

=  -1.7°  X  10-^ 

3s 

=  7.4°  X  10“^ 

%  =  -3°  X  10"'^ 


As  is  shovn  from  the  analysis  of  the  values  of  and  the  shock 

is  very  close  to  a  circular  shock  vave  even  for  large  departure  of  the 
body  from  the  circular  cross  section,  and  the  effect  of  the  terms  corre¬ 
sponding  to  n  =  6,  8,  and  10  is  very  small.  The  velocity  components 
at  the  surface  of  the  body  are  obtained  from  equation  (65)  at  t 

and  the  pressure  distribution  presented  in  figure  13  is  obtained.  In 
the  same  figure,  the  pressure  distributions  obtained  by  using  the  same 
calculated  values  of  ^v^^  ,  (v^^)  ,  and  for  M  =  I.8I  around 

an  elliptical  cone  vith  a  cross  section  having  axes  in  the  ratio  of  1 
to  1.88  are  also  shovn.  The  conical  body  having  an  ellipse  of  axis 
ratio  1  to  3  has  the  same  cross-sectional  area  as  a  circular  cone 
of  =  11^.  Its  pressure  drag  obtained  from  this  calculation 

is  Cj)  =  0.099^  coir5)arison  with  0.12  for  the  circular  cone.  The 

conical  body  having  as  cross  section  an  ellipse  vith  axes  in  the  ratio 
of  1  to  1.88  has  a  drag  coefficient  of  O.IO3,  while  the  equivalent 
circular  cone  of  %  =  10°  30»  has  C-q  -  0.II5.  Therefore,  those 
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calcula'bions  indicaise  tihal:  conical  "bodies  of  circular  cross  section 
have  larger  drag  than  cones  of  elliptical  cross  section. 

The  results  obtained  agree  veil  vith  the  experimental  results, 
also,  if  the  body  shape  chosen  requires  large  values  for  the  ^ 
angle  fields  (v^)  ,  and  any 

other  conical  shape  having  two  planes  of  symmetry  represented  by  an 
equation  (8I)  when  is  10°  can  be  obtained  at  M  =  I.816.  If  the 

flow  has  only  one  plane  of  symmetry,  only  the  terms  in  cos  9,  cos  39, 
and  so  forth,  must  be  considered;  while,  if  no  symmetry  exists,  terms 
in  sin  nS  and  cos  n0  must  be  considered. 

The  flow  fields  defined  by 

obtained  and  given  in  tabulated  form  without  a  large  amount  of  numerical 
work  as  for  circular  cones,  and,  therefore,  the  determination  of  conical 
bodies  can  be  performed  without  difficulty  in  a  very  short  time. 

When  the  shock  shape  is  somewhat  different  from  a  cone  having 
circular  cross  section,  the  basic  flow  field  must  be  different  from  the 
axially  symmetric.  However,  if  the  basic  conical-flow  components  are 
expressed  in  the  form 


Vr  =  Vr  +  Vj.  f(0) 

a  b 


W  =  w^f'(0) 


the  basic  flow  can  still  be  obtained  by  solving  ninnerically  the  equa¬ 
tions  of  motion  in  two  meridian  planes,  and,  therefore,  the  basic  flow 
can  be  determined  exactly.  For  conical  flow  the  linearized  method  can, 
then,  have  wide  application  to  any  form  of  boundary  conditions. 


Flow  Fields  Around  Slender  Bodies  without  Symmetry 

When  the  conical  flow  is  determined,  the  method  of  characteristics 
can  be  applied  to  the  determination  of  slender  bodies.  The  equations 
used  are  similar  to  the  equations  for  circular  bodies  at  angles  of 
attack,  and  can  be  directly  derived  from  those  equations  (reference  2). 

Few  values  of  n  are  required  for  the  determination  of  the  flow 
field,  and  one  set  of  calculations  can  be  used  for  several  bodies  having 
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the  same  hasic  body;  therefore^  the  method  can  he  of  interest  for  prac¬ 
tical  applications. 


Quasi-Two-Dimensional  Flov  Fields 

In  many  general  three-dimensional  flow  fields  of  practical  Interest 
the  flow  is  not  too  different  from  a  two-dimensional  flow^  and,  there¬ 
fore,  the  velocity  field  and  entropy  field  can  he  expressed  as  in  equa¬ 
tions  (1)  and  (2)  with  good  practical  approximation.  Flow  fields  of 
this  kind  are  found,  for  example,  in  wings  having  plan  forms  which  can 
he  considered  close  to  the  two-dimensional  type  with  some  twist  or  a 
variation  of  thickness  distribution  along  the  span.  Flow  fields  of 
this  kind  can  he  considered  also  in  some  problems  in  which  interference 
between  a  wing  and  a  two-dimensional  tall  (downwash  effects)  or  between 
a  two-dimensional  wing  and  a  body  is  considered.  In  all  these  problems 
of  practical  interest  for  the  airplane  design,  the  coxaponent  w  in  the 
direction  of  the  span  of  the  wing  can  be  considered  small;  therefore, 
equation  (1)  can  be  used  and  the  con^onents  u^^  and  v^  depending  on 

the  three-dimensional  effect  can  also  be  considered  small. 

Equation  (17)  expressed  along  the  characteristic  line  of  the  first 
family  =  tan  +  9^^  iii  the  plane  z  =  Constant  is: 


1  dV-i  dcpi  sin^p^  dS-i  Vt  dwn  n  tan  sin  3- 

- i  -  tan  3^  i  +  cp  Ct  =  — i  - ^ 

Vo  dx  o  dx  7R  dx  ^1  1  Vo  1  dz  Vo  cos(cp^  + 

(82) 

while  along  the  second  characteristic  line  >>2  =  “  ^o) 

following  equation  is  valid: 


1  dV-i  dcp-,  sin^3o 

i  ^  +  tan  ^  1 


Vq  dx 


dS- 

dx 


Co  =  ^  i  pQ  Sin  Pq 
•iz  Vo  cos(cpo  -  Po) 


(83) 


where  B-j^,  ^2^ 

streamline  s^. 


C^,  and  are  defined  by  equations  (20).  Along  each 


dwj^ 

dSo  ^z  3z  7EV0 


(81^) 
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Equations  (82),  (83),  and  (84)  permit  the  determination  of  the  flow 
field  hy  relatively  siiiq)le  procedures. 

Consider,  for  exaiiq)le,  a  wing  having  twist,  variable  profile  dis¬ 
tribution,  and  variable  chord,  as  shown  in  figure  l4.  The  wing  can  be 
analyzed  by  means  of  the  linearized  characteristics  method  in  the  fol¬ 
lowing  way:  First,  the  root  and  tip  profiles  are  considered.  Section  a 
and  section  b  have  different  relative  thicknesses  and  chords. 

If  the  variation  from  a  to  b  is  linear,  the  properties  of  a 
two-dimensional  cross  section  at  any  station  c  can  be  obtained  by 
means  of  linear  interpolation  between  the  corresponding  values  at  a 
and  b. 

The  profiles  a  and  b  are  analyzed  by  means  of  two-dimensional- 
flow  theory  and  the  characteristic  net,  and  the  values  of  the  coef¬ 
ficients  B2,  and  C2  are  determined  from  tvo-dimensional 

considerations.  If  entropy  effects  are  neglected  or  incoip)orated  in 
the  linearized  flov^  the  coefficients  B2,  and  C2  can  he 

determined  as  for  the  case  of  two-dimensional  potential  flow  at  each 
point  of  the  axis  (equations  (29),  (30),  and  (3^))  and  are  constant 
along  characteristic  lines  of  the  first  family.  Then  the  linearized 
flow  is  defined  as  the  flow  that  considers  the  three-dimensional  effects 
and  the  entropy  distribution.  Therefore, 


U  =  UQ(x,y,z)  +  aui(x,y,z) 
■V  =  Vo(x,y,z)  +  avi(x,y,z) 
w  =  awi(x,y,z) 


(85) 


where  u^  and  Vq  are  the  potential- flow  solutions  in  the  plane 

z  =  Constant  and  satisfy  the  boundary  conditions  in  the  plane 
z  =  Constant,  U2_  and  are  the  components  due  to  the  presence 

of  w-]_  (and  of  the  variation  of  entropy),  and  a  can  be  a  coefficient, 

for  exaii5)le,  proportional  to  the  twist  distribution  or  to  the  thickness 
variation.  Because  u^  and  Vq  are  functions  of  x,  y,  and  z,  while 

for  the  basic  flow  they  have  been  determined  from  two-dimensional  con¬ 
siderations,  and  are  not  zero;  therefore,  equation  (22) 

oz  oz 


Sw  5w  u  .  ^  V  U  .  V  ft  SV 

- -  S  - - +  - - a:-_n+  —  5  + 

Ssq  hxY  5y  V  V  *  V  Sz 


becomes 


MCA  TN  2515 


47 


or 


SSq  Sz 


+ 


^  _a£_ 

hz  rKVo 


(86) 


where  ^  is  the  variation  of  the  velocity  component  for  the  basic 
dz 

flow. 


When  is  considered  small^  in  all  the  flov  field  the  terms 

^  and  —  can  still  he  neglected  in  the  differential  egua- 

a2  dz  a2  dz 

tions  along  the  characteristic  lines^  and^  therefore,  equations  (82) 
and  (83)  are  still  valid. 


At  each  plane  z  =  Constant,  the  characteristic  net  is  knovn; 
therefore,  the  intersections  of  the  shock  wave  for  the  total  flow  with 
a  plane  y  =  Constant  can  he  determined  from  characteristic  calcu¬ 
lations.  If  0  is  a  point  at  the  leading  edge  of  the  wing  (fig.  14(h)), 
the  shock  wave  at  0  can  he  obtained  from  shock- wave  considerations  and 
from  the  boundary  conditions  because  5  at  0  is  known,  and  at  0  the 
shock  is  two-dimensional.  Therefore,  the  velocity  coi]:5)onents  u^^ 

and  v^  at  0  are  zero,  while  is  given  by  equation  (86).  If  the 

OSq 

plane  x  =  is  assumed  to  he  close  to  the  plane  x  =  x^,  the  charac¬ 

teristic  lines  BA  and  CA  can  he  drawn  in  any  meridian  plane  considered 
for  the  basic  flow. 


At  the  point  0,  is  zero  and,  in  the  nei^horhood  of  0  along 

the  shock  L  the  velocity  can  he  expressed  as 


(87) 


Now,  along  the  shock  wave  the  direction  of  the  velocity  behind  the  shock 
is  related  to  the  intensity  from  the  equations  of  the  shock  wave;  there- 

Scp-| 

fore,  the  direction  cp  =  — ^  along  the  shock  and  the  shape  of  the 

oL 

aVn  SVn 

shock  as  a  function  of  ^  are  determined  ■when  --i  along  the  shock 
is  known. 


1^8 


MCA  TN  2515 


If  the  velocity  Y-^  at  a  point  A  of  the  hody  is 


(88) 


then  the  velocity  at  C  and  B  can 

SV-j 

of  — i  from  equations  (82)  and 
os 


he  determined 
(83),  because 


as  a  function 
the  value  of 


at  A 


is  given  from  equation  (86)  and  is  knovn^  and  the  value  of 
obtained  from  the  value  of  at  A  in  several  planes  z 


6w  - 

—  can  be 
dz 

=  Constant. 


If  in  equations  (82)  and  (83)  the  values  of  and  cp^  at 

B  and  C  are  expressed  by  means  of  equation  (87)  and  of  the  equations 
of  the  shock  waves  which  give  the  coefficients  of  the  expressions 


hx  5L  dx 

6cp-^  q 

5x  6x 

6x  SV-j^  6x 


^  Sv-i 

then  equations  (82)  and  (83)  give  two  relations  between  and  — 

os  oL 

3Vq  ^Vl 

and.  therefore.  — and  — =•  can  be  determined. 

6s  6l 


The  equation  of  the  shock  wave  can  relate  u  as  a  function  of  v, 
or  as  a  function  of  cp^  without  the  necessity  of  the  component  w. 


because  the  component  w  is  proportional  to  the  inclination  q  of  the 
tangent  to  the  shock  with  the  plane  x  =  Constant^  and  and.  cp^  are 


functions  of  cos  q,  but  cos  q 


Ml 


in  the  approxi¬ 


mation  considered  here  (fig.  14(c)).  The  components  u  and  v  at  B 


and  the  position  of  B  having  heen  determined,  in  each  meridian  plane 
z  =  Constant,  the  intersection  of  the  shock  in  the  plane  y  =  Constant 
is  obtained  and  v  at  B  detemined.  Then  a  point  D  is  interpolated 
in  each  meridian  plane  and  the  point  E  is  obtained.  Then  point  F  is 
determined.  In  order  to  obtained  ('W')-ri,  the  streamline  DF*  for  the 

/V  * 


basic  flov  must  be  drawn  and 


Sz 


interpolated  between  F*  and  E. 


In  a  similar  way,  all  the  flow  field  can  be  obtained.  The  line  TT* 
defining  the  plan  form  in  figure  l4(a)  must  be  outside  of  the  Mach 
conoid  from  T.  Because  few  points  along  each  profile  are  required,  the 
largest  amount  of  work  for  such  a  calculation  is  represented  by  the  con¬ 
struction  of  a  basic  characteristic  net  which  permits  obtaining  points 
which  simplify  the  determination  of  w.  By  changing  the  value  of  the 
coefficient  a,  different  thickness  distributions  or  different  twists 
can  be  considered.  The  new  distributions  must  be  obtained  by  changing 
proportionately  the  variation  of  thickness  or  twist  with  respect  to  the 
basic  wing  and  by  varying  in  proportion  the  value  of  a. 


CONCLUDING  EEM/VEKS 


The  method  of  characteristics  for  supersonic  flow  has  been  sln^)!!- 
fied  by  assuming  that  one  of  the  velocity  components  or  the  effect  on 
the  velocity  components  due  to  variation  of  one  physical  parameter  is 
small,  so  that  the  square  of  the  velocity  components  considered  small 
can  be  neglected.  By  means  of  this  simplification,  the  flow  field  can 
be  represented  as  the  superposition  on  a  basic  flow  field  (-vdiich  is  not 
linear  and  must  be  determined  by  the  method  of  characteristics)  of 
linearized  flow  fields  which  are  defined  by  a  differential  equation 
with  variable,  but  known,  coefficients. 

The  calculations  of  these  linearized  flow  fields  can  be  performed 
along  the  characteristic  net  of  the  basic  flow  field.  The  method  has 
been  applied  (a)  to  the  two-dimensional  flow  with  entropy  gradient, 
which  has  been  transfoimed  to  a  basic  potential  flow  on  which  a  line- 
earized  flow  due  to  the  entropy  gradient  is  superposed,  (b)  to  axially 
symmetric  problems  where  conical  or  cylindrical  flows  are  considered 
as  the  basic  flow,  (c)  to  the  determination  of  the  flow  field  around 
cones  or  slender  bodies  without  axial  symmetry,  and  (d)  to  particular 
three-dimensional  flows  which  can  be  simulated  as  a  basic  two- 
dimensional  flow  on  which  three-dimensional  linearized  flows  are  super¬ 
posed.  Application  (b)  permits  obtaining  in  a  simple  way  the  flow 
field  around  bodies  of  revolution  without  using  linearized  theory,  and 
indicates  the  possibility  of  using  tabulated  values  for  such  determi¬ 
nation.  Application  (c)  permits  the  determination  of  flow  fields  not 
yet  determined  by  the  method  of  characteristics.  Any  such  conical 
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flovs  can  be  determined  by  using  tabulated  values  that  can  be  obtained 
as  for  cones  of  circular  cross  section  at  small  angles  of  attack.  The 
application  in  (d)  can  be  of  interest  for  vings  of  approximately  two- 
dimensional  form  having  twist  or  thickness  variation  along  the  span  and 
to  interference  problems. 
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Fig\u:e  4.-  The  basic  net  when  conical  flow  is  assumed  as  the  basic  flow 
,  for  the  linearized  characteristics  method. 


Characteristics  method 
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Figure  7.-  Comparison  between  results  of  the  calculation  by  the  character¬ 
istics  method  and  by  the  linearized  characteristics  method  obtained 
at  M  =  3.016. 


igure  9*“  Determination  of  linearized  flov  of  constant  intensity  superposed 

on  conical  flov. 


10.-  Application  of  the  linearized  characteristics  method  having 
constant  disturbance  superposed  on  a  conical  basic  flow. 
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